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We derive an integral convex combination of product states for a range of separable Werner 
states. Our method consists of expanding the sought-after local density operators in terms of 
Wigner operators. For dimension d = 2, our decomposition holds for the whole separable range of 
Werner states, while for d > 2 it is valid for a subset of separable Werner states. We illustrate the 
general method with the explicit examples d = 2 and d = 3. 



I. INTRODUCTION 
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5^ ' Composite quantum systems are among the foci of quantum information theory. Bipartite quantum states, i.e. 
states of quantum systems that consist of two subsystems A and B, can be classified according to their property of 
being entangled or separable. In 1989, R. Werner proposed a physically meaningful definition of separability of a 
quantum state [1], namely 
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PtPA ® Pb i (!) 



where are probabilities, i.e. Pi > and YliPi = 1. In the integral version of this convex combination of product 
states, the discrete index i is replaced by some continuous variable A, such that the integral product decomposition 
■ reads p sep = J dXp(X)pA(X) <8> Pb(X), where p(A) > and J dXp(X) = 1. Given a separable state, its product 
£C) I decomposition is not unique, and in general, it is a hard task to find such a decomposition. If one does not know 
whether a given state is separable or entangled, it is difficult to prove whether a separable decomposition does or does 
not exist. This is the origin of the entanglement vs. separability problem [2, 3]. 

In his seminal paper [1], Werner also introduced a certain family of states that is nowadays referred to as Werner 
states. Due to their symmetry properties, they play an important role in several contexts, e.g. in entanglement 
Oh! purification [4], entanglement properties for states with white noise, and the possible existence of bound entangled 
t_L ■ states with non-positive partial transpose [5] . 

For qubits, a separable decomposition of Werner states has been found in [6]. A different decomposition was 
recently derived in [7], again for 2-dimcnsional subsystems. Here, we find a separable decomposition of Werner states 
in arbitrary dimensions. For dimension d = 2, our decomposition is different from both decompositions in [6] and [7]. 



II. WERNER STATES 



Werner states [1], in the following denoted as pw, are a class of mixed states for bipartite quantum systems (where 
each of the two subsystems has dimension d), which are invariant under the transformations U ®U , for any unitary 
U, i.e. pw = (f<£> U)pw{U' (g> W). The family of Werner states, characterized by one parameter /, is (in the original 
notation) given by the density operator 

Pw = j^2((d-f)t + (df-l)V) , (2) 

which acts on the d x d-dimensional Hilbert space Ha <£> Hb that is spanned by the state vectors of subsystems 
A and B. Here, 1 is the identity operator on Ha <8> Hb, and by V we denote the swap operator, which acts as 
V| 0)^4 <8> \ iP)b = \^)a ®\4>)b- Positivity of pw implies for the parameter / = tr(pw'V) that — 1 < / < 1. The 
Werner state pw is separable, i.e. classically correlated, iff < / < 1. In the following, we will regard the Hilbert 
space of a subsystem as the state space of a particle with spin j, where d = 2j + 1. The basis states are denoted by 
\jm), where m = j, j — 1, ... — j. It is known [8] that V has the following form in terms of the Wigner operators 
[8]: 

V = EE (-l) q T^Tl ql (3) 

k=0q=-k 
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where denotes the (2j + 1) x (2j + 1) matrix with the elements given by [8] 



I 2k ■ 1 y 2 ,, , 



(jmi|T g fc |jm 2 ) = 



2j + 1 



(4) 



The coefficients C^^ m2 are Clebsch-Gordan coefficients. We use the notation | JMjij 2 ) — J2 mi m 2 ^mlm 2 u\ h m i) ® 
IJ2W2), where J is the total angular momentum, M its third component, ji is the angular momentum of particle i, 
and mi its third component, with mi + rri2 = M. Note that the (2j + l) 2 operators T£ generate the group U (2j + 1). 
The matrices obey the orthogonality relations 



1 <I2 



tr (T^.(T q k 2 ^)=5 klk2 6 qiq 
and are traceless for k ^ [8] . For k — we have 

t ° = (2.? + f r i/2 1. 

Using Eq. (3), we can rewrite the Werner state (2) in terms of the Wigner operators: 



Pw 



Y k=0q=-k J 



(5) 



(6) 



(7) 



where 



III. DECOMPOSITION OF SEPARABLE WERNER STATES 

In general, a separable state can be written in the integral decomposition [1] 

Psep = J dXp (A) [ PA (A) ® p B (A)] , (8) 

p(A)>0, [d\p(\) = l, 

and (A) > 0, /?b (A) > represent density operators of subsystems A and -B. Here, A has to be understood as a 
symbol for one or more continuous variables. 

We can express an arbitrary density operator p (A) of a spin-j particle via the operator decomposition 

2j k 



where the expansion coefficients y kq are given by 

y kq (X)=tr{p(\)T^) , 

due to the orthogonality relation (5). 

Using this decomposition for both Pa(X) and pb(X) we can rewrite a separable Werner state (7) as 



(9) 



Pw 



= J d\p(X) 
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(10) 



A comparison of (7) with (10) leads to the following condition on y kiqi (A) and y k2q2 (A): 

J d\p(\) y kiqi (A) y k2q2 (A) - S kltk2 5 qi _ 92 . 



(11) 
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We can conclude that y kiqi (A) with i = 1, 2 are orthogonal, with a weight function p (A). Remember that the Werner 
state on the left-hand side of Eq. (10) depends on the parameter /, which is not explicitly written here. Thus, the 
expansion coefficients yk iqi (A) and y k2q2 W an( i the weight function p (A) will in general be functions of / as well. 

It is possible to satisfy the orthogonality condition in Eq. (11) in the following way: let us consider A as a symbol 
for two parameters. By defining A = {9, ip} with < 9 < n and < <p < 2ir and p (9, <p) — 1/A-k [7] the relation (11) 
takes the form 

j d9 J d(psm9y kiqi (9,ip)y k2q2 (6,ip) ~ S kl ,k 2 5 gi ,- q2 . (12) 
o o 

Thus we can use for y kq (9, <p) the spherical harmonics Y q (9, <p). One can readily check that the decomposition (10) 
can be rewritten as 

Pw^-^Jdojd^ sin 9 \p A (9, y) ® p B (0, ip)] , (13) 
o o 

where 

« <9,<p ' = W+T 1 + wTU%T^ f&*rc ( " * <14) 



j 2j k 

PB {6, <p) = —— 1 + E E ^ q Y q k {9, <p) T q k , (15) 

■> k=lq=-k 

with r]k q being arbitrary parameters. We note that hermicity of pa.b {0, <p) implies that rj^ q = rjk. - q - We have used 
the properties of the spherical harmonics Y q (9, ip) , namely Eqs.(12) and Y q * (9, <p) = (— l) q Y* q (9, ip). The above 
combination of product states for pw still contains the freedom in choosing the parameters r]k q - Note, however, that 
at this point we have not yet shown positivity of the local operators, and the notation pa and ps is merely suggestive. 
The reader will have noticed that the prefactors of the sums in expressions (14) and (15) are not identical. We have 
chosen this asymmetric decomposition on purpose, as will be explained below. The decomposition (13) could have 
equally well been formulated in a symmetric way by defining 



k— lq— — k 



I ■ 2 j k 

B V ' 2. 7 + l y (2. 7 + 1) (2.7 + l) 2 - 1 ^™ fc 

In this version p^ {9, ip) and p$ {9, ip) have a nice symmetric form, however their positivity conditions are quite 
difficult to analyze. It is obvious that the separability property of the Werner state does not depend on the form 
of the local density operators, therefore we choose to use the local operators (14) and (15), because they allow to 
determine the positivity constraints in an easier way. 

If all eigenvalues of the local operators (14) and (15) are positive, we found a valid separable decomposition of the 
Werner state pw in Eq. (13). The calculation of the eigenvalues of the local density operators for j > 1 is a difficult 
task. We now simplify the analysis by choosing r) kq = r\ k . Since 

Y q k (6,<p) = ] l™^D k qQ *(e,<p) , (18) 

where D qm {9, ip) is the Wigncr rotation matrix, defined via the transformation of T q as a tensor [8], i.e. 

k 

U : T k m -+ UTW = E Dk qm V) T q k ■ (19) 

q— — k 
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Inserting Eq. (18) into Eqs. (14) and (15), using Eq. (19) and the fact that (8, ip) = D k qm (6, —tp), we can rewrite 
the local density matrices in the form 



2j 



(20) 



PB {6, tp) = U B {6, - v ) (^J-fJ 1 + JjlkV^k + lT^j U B {6, -ip) . (21) 

Here, U (6, ip) denotes the unitary irreducible representation of the SO (3) group on the state space spanned by \jm), 
with m = j, j — 1, —j, and is defined as 

U (9, tp) = exp (—itpj z ) ■ exp (—i9J y ) . 

Note that in the unitary operator we have already omitted the third Euler angle 7, by dropping the factor exp (17 J z ), 
because we are using the eigenbasis of J z . 

Thus, we can rewrite the Werner state (2) in the form 



-rr 2-rr 



Pw = ^Jd6 J dtp sin e [U A (0, -<p) ® U B (0, -<p)] [pA ® P b] [Ua {0, -<p) f ® U B (0, -<p) f } , (22) 




where 

2j 



1 V , 

PB = ^—1 + 5>V2fcTlT fc . (24) 

It is easy to see that P a and ps are diagonal matrices. To show this one uses (4) and the properties of the Clebsch- 
Gordan coefficients, namely C^ qm2 = if mi + g ^ m 2 . Thus the diagonal elements, i.e. eigenvalues of P a and /9b, 
are given by 

A - <"> - wrt + wW^rl* ' (24 + 1)cS " (25) 

and 

A ro (p B ) = ^ + ^=^( 2fc + 1 ) C ^' m = j,i-l,...-j. (26) 

If all A TO (pa) and A TO (ps) are positive, we found a separable form of the Werner state (2). We note that positivity 

j 

of A m (p A ) and A m (p B ) for all m implies that ^ A m (pa) A m (ps) > 0. Owing to the fact that [8] 

m=-j 

E ^ol^ol = and E ^ = for fc ^ 0, (27) 

m——j m——j 

3 

one finds ^ A m (p^) A m (p B ) = /• Thus / > is a necessary condition for positivity of p A and ps. It remains to 

m=-j 

show the sufficient conditions for positivity of p^ and ps- At this point it becomes clear why the asymmetric form, 



■5 



chosen above for pa and pb, is advantageous: we can determine the free parameters r/k from the positivity condition 
for pB and then find the range of / for which pa is positive. 

We note that for the simplest case of spin j = 1/2 the present decomposition is, due to the isomorphism between 
SO(3) and SU(2), equivalent to 

p w = JdU [U®U]- [ PA ® pb] ■ [U ® U] f , (28) 

where the integral is extended to all unitary operators acting on the two-dimensional Hilbert space, with J dll = 

1 and dll representing the standard Haar measure on the group SU(2). Hence, we expect that the inequalities 
A±i/2 {pa.b) > will yield < / < 1. For higher dimensions, however, there is no such simple argument, and it 
turns out that for higher spins the decomposition (22) is a separable decomposition, i.e. pa b > 0, only in the range 
< / < f (j) , where f (j) < 1. 

The presented method can also be used to find a separable forms for more general states e.g. states which are 
invariant under product representations of the group SO (3) of three-dimensional rotations, see [9]. 



IV. EXAMPLES 



Let us consider explicitly the two lowest dimensions, namely d = 2 (spin j = 1/2) and d — 3 (spin j = 1). 



A. Casej = l/2 

Here, the index k in Eqs. (25) and (26) takes one value, k = 1, and we have one free parameter 771 . Thus, for qubits 
we have the following inequalities: 



A1/2 (pb) = l( 1 + vfon) > 0, 

A-1/2 (Pb) = \ (l - V6771) > 0. 



This leads to the condition 



Positivity of X±i/2 (pa) translates to the constraints 



The solution of the inequalities (29), (30), and (31) is 

\m\ < -7= and < / < 1. 
V6 

Thus, as mentioned above, for j = 1/2 our decomposition is valid for the whole separable range of the Werner state 
family. 



B. Case j = l 

For j = 1, we have k — 1,2. We can use two free parameters 771 , 772 . Positivity of A m (pb) for m = — 1,0, +1 implies 
the following inequalities: 
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Ai (pb) = l( 2 + + V30m) > , (32) 

A_i (ps) = - (2 - 3V6?7i + V30r; 2 ) > , 
Ao (pb) = l( 1 ~ > . 



The eigenvalues of X m (pa) read 



w , 1 ( Vi (V30 (3/ - 1) + 48>7 2 ) + 3V6 (3/ -l)r, 2 \ 



J_ / Vi (V30(3/ - 1) + 48^) - 3V6 (3/ - l)m 
144 I 772^71 



A-i (pa) 

After a lengthy calculation we conclude that positivity of pa holds for 

o</<|, 

i.e. for the case j = 1 our decomposition is valid only within a certain range of the parameter /, and not for the 
whole separable interval < / < 1. Note that our separable decomposition holds for Werner states that are "close" 
to entangled states. 

Fig.l shows the solutions of inequalities A r „ (pa.b) > for different values of / , for spin j = 1. One can see that 
the common area of filled regions vanishes for / > | . 



V. CONCLUSIONS 



In summary, we have used an angular momentum approach to find in principle a separable decomposition of Werner 
states in any dimension. Our main idea is to express the local density operators and the global operator in terms 
of Wigner operators. The eigenvalues of the local operators are found to be certain combinations of Clebsch-Gordan 
coefficients. Some free parameters allow to guarantee positivity of the local density operators. For dimension d = 2 
our decomposition holds for all separable Werner states, while for higher dimensions it is only valid for a certain range 
of separable Werner states. The size of this range depends on the dimension. We verified that our decomposition is 
valid for / = up to dimension d = 5. It is still an open task to find a product decomposition of all separable Werner 
states in any dimension. 

We would like to thank R. Werner, M. Kleinmann and T. Meyer for informative discussions. 
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FIG. 1: Case spin j — 1: Solutions of inequalities (32) and (33) (filled regions) for different values of / = 0, |, | and 1. 



